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Abstract 

We use an accurate implementation of density functional theory (DFT) to calculate the zero- 
temperature generalized phase diagram of the 4d series of transition metals from Y to Pd as a 
function of pressure P and atomic number Z . The implementation used is full-potential linearized 
augmented plane waves (FP-LAPW), and we employ the exchange-correlation functional recently 
developed by Wu and Cohen. For each element, we obtain the ground-state energy for several 
crystal structures over a range of volumes, the energy being converged with respect to all technical 
parameters to within ~ 1 meV/atom. The calculated transition pressures for all the elements and 
all transitions we have found are compared with experiment wherever possible, and we discuss 
the origin of the significant discrepancies. Agreement with experiment for the zero-temperature 
equation of state is generally excellent. The generalized phase diagram of the Ad series shows 
that the major boundaries slope towards lower Z with increasing P for the early elements, as 
expected from the pressure induced transfer of electrons from sp states to d states, but are almost 
independent of P for the later elements. Our results for Mo indicate a transition from bcc to fee, 
rather than the bec-hep transition expected from sp-d transfer. 
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I. INTRODUCTION 



The transition metals are among the most intensively studied families of elements, and 
there has been an effort going back several decades to map and interpret systematic trends in 
their properties. Some of these trends are of fundamental importance to our understanding of 
the energetics and electronic structure of transition metals. Examples include: the parabolic 
dependence of cohesive energy and bulk modulus on d-band filling; 1 the sequence of most 
stable crystal structures associated with increasing band filling at ambient pressure;- and 
the pressure induced increase of <i-band width going roughly as the fifth inverse power of 
the atomic volume.- However, our knowledge of transition-metal systematics is still very far 
from complete, as is evident from the current major controversies over high-pressure melting 
curves.-^ 1 ^ Even at low temperatures, there are sizable gaps in the map of transition-metal 
phase diagrams, and new crystal structures continue to be discovered.- We report here 
a systematic investigation of the zero-temperature phase diagram of all the Ad transition 
metals over a wide range of pressures, based on one of the most accurate implementations 
of density functional theory (DFT) currently available. 

There have, of course, been very many previous detailed studies of transition metals based 
on DFT, some of which investigated the relative stability of different crystal structures 
at high pressures. However, most previous work has been designed to address particular 
questions relating to particular metals. Here, by constrast, the aim is to obtain a coherent 
overall view of an entire transition-metals series. Specifically, using Z to denote atomic 
number, we wish to use DFT to map out the generalized (P, Z) phase diagram of the Ad 
series at T = K at pressures P up to ~ 500 GPa£. In order to substantiate the accuracy 
of the calculations, we shall compare our calculated results for P as a function of volume V 
with all available experimental data. The calculations are all performed using the FP-LAPW 
implementation of DFT (full-potential linearized augmented plane waves) which for 

a given exchange- correlation functional is among the most accurate ways of calculating the 
total energy. We use the exchange-correlation functional recently developed by Wu and 
Cohen,— which appears to reproduce the experimental energetics of Ad transition metals 
more accurately than other functionals.— 

There are several motivations for wishing to obtain the (P, Z) phase diagram of a 
transition-metal series at T = K. One motivation is the possibility of discovering hitherto 
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unknown crystallographic phase boundaries for particular metals. Another motivation is to 
probe our basic understanding of transition-metal energetics. At ambient pressure, the most 
stable crystal structures of the Ad and 5d series follow the sequence hexagonal close packed 
(hep) to body-centred cubic (bec) to hep to face-centred cubic (fee) with increasing Z, and 
the 3d series nearly follows this sequence, the deviations being caused by magnetic effects.-- 
This sequence is entirely explained by band energies, and can be reasonably well reproduced 
using a canonical ci-band tight-binding model with some modifications due to hybridization 
with sp bands.— The main effect of pressure is to shift the relative energies of the (i-band 
centroid and the bottom of the sp band, and hence increase the <i-band filling.— This leads 
one to expect that the same sequence of structures will be found at high P, but with the 
boundaries shifted to the left (towards lower Z) .— »21 We shall confirm that this is the case, 
though the reality is slightly more complicated. 

In addition to the motivations we have mentioned, there is another that is important to 
us. Current controversies over the high-pressure melting curves of transition metals stem 
from apparent disagreements between melting temperatures deduced from static compression 
and shock experiments, these disagreements amounting to several thousand K at megabar 
pressures.- 1 ^ DFT calculations of transition-metal melting curves^ 1 ^ support the correctness 
of the shock results, and one of the proposed resolutions of the conflicts is that the transitions 
identified as melting in static compression are in fact solid-solid transitions. It would clearly 
help to reduce the confusion if we had a better understanding of the generalized (P, T, Z) 
phase diagram of transition- metal series 9 . At present, we have a fairly complete knowledge 
of this only at low P, and, to a lesser extent, at low T. We see the present attempt to 
complete the low-T (P, Z) diagram as an essential step towards mapping the full (P, T, Z) 
diagram. 

Earlier DFT work has already given much information about the energetics of transition 
metals at T = K. Particularly important here is the work of Pettifor— i22j2£ and others, 
which provided a systematic understanding of the structural trends across all the transition- 
metal series in terms of the electronic densities of states of different crystal structures. Their 
work also gave important insights into the pressure induced transfer of electrons from sp- 
states to d-states. Also important was the first-principles work of Skriver— on structural 
trends across transition metal series at zero pressure. The pioneering work of Moriarty, 
Johansson and others on the high-pressure energetics of transition metals will be cited below. 
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The information from all this previous work could be assembled to produce a substantial 
part of the (P, Z) diagram, but there would be inconsistencies from the use of different DFT 
implementations and exchange-correlation functionals. Here, we avoid such inconsistencies 
by using a single high-accuracy method for all the calculations. 

The remainder of the paper is organized as follows. In Sec. [TT1 we summarize the essential 
ideas underlying the methods used in this work and the results of our convergence tests. Then 
(Sec. IHIj) . we present for each of the transition metals from Y to Pd our calculated energy 
differences between different crystal structures and the equation of state (EOS) P(V) for 
pressures from ambient to ~ 500 GPa, comparing our results with experiments and previous 
calculations where possible. At the end of Sec. 1111} we present the zero-temperature (P, Z) 
phase diagram of the 4d transition-metal series. Discussion and conclusions are in Sec. IIVI 

II. TECHNIQUES AND TESTS 

There are several implementations of DFT that can be used to calculate the total energy 
per atom of a crystal, including pseudopotential techniques,— the projector augmented wave 
technique,— 1 ^ the full-potential augmented plane-wave (FP-LAPW) technique,— etc. 
We have chosen to use FP-LAPW here, because for a given exchange-correlation functional 
it can give values for the total energy that are closer to the exact value than most other 
implementations. There are a number of technical parameters in FP-LAPW that control 
convergence towards the exact value. We recall in this Section the main ideas of FP-LAPW, 
summarize the parameters that control convergence, and report tests that guide our setting 
of these parameters. The exchange-correlation functional used for all the present work is 
that due to Wu and Cohen.— Later in this Section, we outline briefly what this functional is, 
and we report tests indicating that it should be particularly accurate for present purposes. 

A. Full- potential linearized augmented plane waves 

In APW methods, space is divided into non-overlapping muffin-tin spheres (radius -Rmt) 
centred on the atoms, and the interstitial region between the spheres. Each Kohn-Sham 
orbital is represented as a sum of radial functions multiplied by spherical harmonics up to a 
maximum angular momentum Z™^ x in each sphere, and as a sum of plane waves (maximum 
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wavevector K max ) in the interstitial region; the coefficients are adjusted to achieve continuity 
at the sphere boundaries. The Kohn-Sham eigenstates are divided into core and valence 
states; core states are treated as being non-zero only inside the muffin-tin spheres, while 
valence states extend over all space. The radial functions used as basis sets are solutions of 
the Schrodinger equation inside the atomic spheres, and are therefore energy dependent. In 
linearized APW (LAPW), the energy dependence is treated in a linear approximation. In 
full-potential implementations (FP-LAPW), the Kohn-Sham potential is also represented as 
a sum over angular momentum functions within the spheres (maximum angular momentum 
^poT) an d a sum over plane waves (maximum wavevector G ma x) in the interstitial region. As it 
stands, this scheme is not accurate if there are semi-core states, i.e. low-lying states that are 
not adequately treated by linearization. One solution to this problem is to linearize the basis 
functions using different reference energies for valence and semi-core states. Alternatively 
one can augment the basis set with "local orbitals", this procedure being known as the 
"APW+lo" method. We use the full-potential version of the latter method here in order 
to treat all the valence and semi-core states. A fully relativistic treatment is used for core 
states, and a scalar relativistic treatment for all other states. The wavevectors k for which 
the Kohn-Sham equations are self-consistently solved must be sampled over the Brillouin 
zone of the crystal, and the well-known Monkhorst-Pack fc-point sampling scheme is used 
for this purpose. 29 All the calculations were performed using the WIEN2k code.— 

Throughout the present calculations, we aim to achieve convergence of the total energy 
to within 1 meV/atom with respect to all the parameters we have just mentioned. All states 
having principal quantum number n < 3 are treated as core states, and all higher states as 
valence states. This choice is based on the fact that the states with n = 3 lie at least 200 eV 
below the Fermi level so they will not respond significantly to compression even at megabar 
pressures. Provided the muffin-tin spheres do not overlap, and provided the calculations 
are fully converged with respect to all parameters, the choice of -Rmt should not affect the 
results, but it does affect the efficiency of the calculations. We used the default setting of 
-Rmt provided by the WIEN2k algorithm, which ensures that -Rmt varies appropriately as the 
volume per atom is varied. To determine the settings required for the technical parameters 
-K'max) C max , /™^ x an d ^o°b> an d f° r ^-point sampling, we have conducted systematic tests on 
bcc, fee and hep Mo; we assume that the settings that give the required accuracy for Mo 
will also serve for all the 4d transition metals in all the structures of interest. 
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Our tests on Mo indicate that the choices RMTK max = 9.5, G max = 18 bohr -1 , /™^ x = 10 
and /p^ x = 4, together with 14 x 14 x 14 fc-point sampling (14 x 14 x 7 for hep), give 
convergence to within ~ 1 meV/atom. To demonstrate that this is the case, we have set all 
of these parameters except one to the values just quoted, and studied the effect of varying 
the free parameter. The values of total energy per atom E tot obtained when we vary in turn 
-RMT-Kmax, G max , , l 1 ^ and fc-point sampling are reported in Table I , from which we 
see that the required convergence is indeed achieved with the quoted values. These are the 
values used for all the calculations reported later. We note that, although we have paid 
careful attention to convergence, there remain small errors due to the linearization inherent 
in the FP-LAPW scheme. 

B. Wu-Cohen exchange-correlation functional 

The Wu-Cohen exchange-correlation functional is a particular form of generalized gradi- 
ent approximation (GGA). In GGAs, the total exchange- correlation energy is the integral 
over the volume of the system of a density of exchange-correlation energy, this density be- 
ing expressed in terms of the number density n(r) of electrons as n(r)e xc (n(r), s(r)). Here, 
e xc , the exchange-correlation energy per electron at position r, depends not only on n(r) 
itself but also on the magnitude of its dimensionless gradient s(r) = | Vn|/[2(37r 2 ) 1//3 n 4,/3 ]. 
The energy e xc is expressed as e xc (n(r), s(r)) = e™ lf (n(r))F xc (s(r)), where e^(ri) is the 
exchange-correlation energy per electron in the uniform electron gas of density n, and F xc 
is the so-called enhancement factor. 

In the widely used Perdew-Burke-Ernzerhof (PBE) form of GGA,— a parameter-free 
formula expressing the dependence on s of F xc (s) was derived by requiring that certain 
exact conditions be satisfied. The WC functional takes the same GGA form as PBE for the 
correlation part of F xc , but modifies the exchange part. The modification was done with the 
intention of improving the functional for condensed matter, at the expense of a somewhat 
worse description of atoms and molecules. Specifically, the modification was based on the 
idea that the exchange part of F xc should be constructed so as to reproduce the fourth-order 
gradient expansion of the exact exchange-energy functional of the electron gas in the limit 
of slowly varying density — The detailed form of the exchange enhancement factor in the 
WC functional can be found in the original paper, 14 which also presents the results of tests 
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TABLE I: Tests on convergence of FP-LAPW total energy per atom E to t (eV units) for bcc Mo. 
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suggesting that the functional gives better results than LDA and PBE for the equilibrium 
lattice parameter and bulk modulus of a range of materials. Subsequently, tests of WC were 
reported on a much more extensive set of materials using the FP-LAPW implementation of 
DFT. 15 These tests supported the claim that WC is on average more accurate than LDA or 
PBE, but they showed that the improvement over these other functionals is by no means 
uniform. However, for Ad transition metals at ambient pressure, the WC predictions of 
equilibrium lattice parameter and bulk modulus are, with a few exceptions, better than 
both PBE and LDA. 

To confirm that the quality of WC is maintained at high pressures, we compare in Fig. [1] 
our calculated results for pressure P as a function of volume V at zero temperature for 
Mo obtained with WC, LDA (Ceperley-Alder)^ and PBE with experimental results up to 
P = 300 GPa. The experimental measurements were performed at room temperature but 
were corrected for thermal effects so that they refer to T = K. 35 The comparisons show 
that WC is in almost perfect agreement with experiment over the entire pressure range, 
while both LDA and PBE show significant deviations. Comparisons of calculated P(V) 
curves with experimental data for other Ad transition metals presented in the next section 
further confirm the quality of WC. 



III. RESULTS 



For each of the Ad transition metals from Y to Pd, we have calculated the total energy 
at a closely spaced series of atomic volumes for the body-centred cubic (bcc), face-centred 
cubic (fee) and hexagonal close-packed (hep) structures. In some cases, we have also studied 
other crystal structures. The pressure as a function of volume is then obtained for each 
structure by fitting the total energy results with a 3rd-order Birch-Murnaghan equation of 
state.™ From these results, we straightforwardly obtain the most stable crystal structure 
at each pressure, as well as the pressures of the structural transitions and the volumes 
of the coexisting phases at these transitions. In the following, we present first for each 
element the difference between the energy per atom in each structure and the energy in 
bcc structure, as a function of V. We then present P(V) of the stable structures over the 
pressure range from ambient to typically 500 GPa. Where possible, we compare P(V) and 
the transition pressures with experimental measurements, and we indicate the relation with 
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FIG. 1: Zero-temperature equation of state of Mo up to P ~ 400 GPa calculated with GGA(WC), 
GGA(PBE) and LDA(CA) exchange-correlation functionals. Dots show experimental data^ 
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TABLE II: Calculated volumes V±, V 2 (units of A 3 /atom) and transition pressures Pt (units of 
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FIG. 2: Top: Energy differences AE of the hep, aSm (C19), dhep (A3') and fee structures with 
respect to bee in Y as function of volume per atom. Bottom: Equation of state of Y up to 
P ~ 700 GPa. 



A. Yttrium 



Experimentally, Y has the hep structure at ambient pressure, and with increasing pressure 
passes successively to the a-samarium (aSm), double hep (dhep) and fee structures.— 
This is the same sequence of structures observed in the lanthanides under pressure, as 
has often been discussed. 40 Earlier DFT calculations on Y have correctly predicted these 
structures,— and also indicated a transition to bee at very high pressure.— Fig. [2] reports 
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our calculated energies of the hep, aSm, dhep and fee structures relative to bec as a function 
of V. This shows the sequence of stable structures in the experimentally observed order. 
The resulting EOS is reported in Fig. [21 and the transition pressures and volumes are given 
in Table II . Although the calculated sequence of crystal structures is the same as that 
observed experimentally, we significantly underestimate the transition pressures by at least 
10 GPa. The difference is considerably greater for the transition aSm — > dhep, but it is clear 
from Fig. [2] that the energy curves of these two structures follow each other so closely that 
the prediction of this transition pressure is likely to be very challenging. Possible reasons 
for our underestimation of the transition pressures are discussed in Sec. HV] . 

We also report in Table HT1 the results of very recent DFT work on Y up to ~ 150 GPa 
by Lei et a/. 41 , as well as those of the older work of Melsen et alM. The hcp-aSm and 
aSm-dhcp transition pressures Pt of Lei et ai, who used FP-LAPW based on LDA, are 
within 2 GPa of ours. It is known that some tranverse phonons of fee Y are unstable at 
pressures below 90 GPa.— Based on this fact, Lei et ai. analyze a distorted fee structure 
which they find to be more stable than fee. In order to save effort, however, we consider 
here the fee structure directly. It is striking that the fcc-bcc Pt of Melsen et ai, obtained 
from FP-LMTO calculations based on LDA is below ours by over 200 GPa. It appears to 
us that their treatment of semi-core states was much more approximate than the one used 
here, and this may account for the large difference. The very small volume change of only 
~ 0.04 A 3 /atom that we find for this transition will magnify the effect of any technical errors 
in the prediction of P t . We return to this question in Sec. HVl 

B. Zirconium 

Experimentally, the stable structure of Zr at ambient pressure is hep, but measurements 
show a transition to the u phase at the rather low pressure of 2.8 GPa, followed by a 
transition to bee at P ~ 30 GPa.— Our calculated energies of the u, hep and fee structures 
relative to bee as a function of volume (Fig. [3]) show the same sequence, and our calculated 
P(V) in the u structure is in close agreement with experimental data. However, as in the case 
of Y, our calculations appear to underestimate the transition pressures (see Table III). We 
find that at P = 0, T = 0, the u structure is most stable, the transition between hep and u 
occurring at the negative pressure of —8.4 GPa. Our calculated transition pressure between 
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TABLE III: Calculated volumes V\, V% (this work, units of A 3 /atom) and calculated and experimen- 
tal transition pressures Pt (units of GPa) for structural transitions in Zr. Calculated Pt values are 
listed according to the exchange-correlation function (WC, LDA or PBE) used. Sources: a: present 
work; b: Ostanin et aL 46 ; c: Schnell and Albers 4 ^; d: Jona and Marcus 4 ^; e: Zhao et al. 44 ; f: Hui et 
alM. 

lu and bcc of 22.2 GPa is also significantly below the experimental value of 30(2) GPa.— 
The underestimation by ~ 10 GPa is similar to what we found for Y, and will be discussed 
in Sec. ED 

Previous theoretical estimates of the u;-bcc transition pressure P t show considerable vari- 
ation (Table [TTTT) . It seems that Pt from LDA is lower than that from PBE by ~ 10 GPa, 
but there are significant differences between results obtained with the same DFT functional. 
The results of Jona and Marcus^ seem to be seriously out of line with both experiment and 
other calculations. We shall comment further in Sec. IIVI 

C. Niobium 

Experimentally, the observed structure of Nb is bcc over the entire pressure range from 
ambient to ~ 145 GPa. Our total-energy results (Fig. HJ) show that hep, fee and u are all 
less stable than bcc by at least 0.2 eV/atom over the pressure range up to over 400 GPa. 
A possible phase transition at pressures above 400 GPa may be suggested by the down- 
turn in the energy difference between hep and bcc at ~ 400 GPa, but we have not pursued 
this possibility. The agreement of our calculated P(V) curve with the recent experimental 
data of Kenichi and Singh™ obtained at room temperature is excellent. Our calculated 
equilibrium volume Vq = 17.94 A 3 /atom, agrees closely with the experimental value at 
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FIG. 3: Top: Energy differences AE of the hep, uo (C32) and fee structures with respect to bec 
in Zr as function of volume. Bottom: Equation of state of Zr up to P ~ 600 GPa; dots represent 
experimental data from Ref. 44] . 



ambient conditions^ of 17.98 A 3 /atom. 



D. Molybdenum 

The properties of Mo have been intensively studied both experimentally and theoretically 
over the last decade. We mentioned in the Introduction the continuing controversy over its 
high-pressure melting curve. At low temperatures, experiments show that Mo has the bec 
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FIG. 4: Top: Energy differences AE of the oj, hep and fee phases with respect to bcc in Nb 
function of volume. Bottom: Calculated EOS of Nb compared with experimental data of Ref. 
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structure for all pressures up to ~ 400 GPa. There have been at least seven previous DFT 
studies on its phase transitions at higher pressures, and all agree that there is a transition 
to either hep or fee, but there is no consensus about which high-P form is more stable. 
There is also considerable variation of the predicted transition pressures. Our calculated 
energy differences fee - bcc and hep - bcc (Fig. |5]) show clearly that bcc is the stable phase 
up to a pressure of nearly 660 GPa. Above this, both fee and hep fall below bcc, but fee 
is below hep, and the fee - hep difference becomes more negative with increasing pressure. 
Our calculated P(V) curve is in excellent agreement with experimental data over the whole 
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FIG. 5: Top: Energy differences AE of the fee and hep phases with respect to bee in Mo as function 
of volume. Bottom: Equation of state of Mo compared with experimental data.— 

pressure range up to 300 GPa for which there are data.— 

We summarize in Table [IV] all DFT calculations on the high-P transition in Mo, noting 
the predicted high-P structure and transition pressure P t to that structure, together with 
the DFT implementation and exchange-correlation functional. In general, we believe that 
work in which careful attention has not been give to convergence with respect to all technical 
parameters and the treatment of semi-core states must be regarded as less reliable. However, 
it remains difficult to draw conclusions about the causes of the very large differences be- 
tween predictions, in some cases based on exactly the same exchange-correlation functional. 
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GGA(WC) 


660 



TABLE IV: DFT predictions of the high-pressure phase of Mo, with DFT method and exchange- 
correlation functional used, and transition pressure P± (GPa units) from works cited in first column. 
Forms of LDA are due to Hedin-Lundqvist 57 (HL), Perdew-Zunger 58 (PZ) and von Barth-Hedin 59 
(BH), and of GGA to Perdew-Wang22 (PW) and Wu-Cohen 14 (WC); the form of GGA used by 
Jona and Marcus was unspecified. 

Further comments will be made in Sec. UVl 
E. Technetium 

Tc has the hep structure at ambient pressure, but little seems to be experimentally known 
about its behaviour under pressure. Our results for the energies of the fee, hep and simple 
hexagonal structures relative to bec (Fig. [6]) show that hep is the most stable up to pressures 
of at least 500 GPa, and there is no indication that any of the other structures will become 
more stable at higher pressures than this. The calculated P(V) curve of hep Tc is reported 
in Fig. El but there appear to be no experimental data to compare it with. Our value 
for the equilibrium volume of 14.21 A 3 /atom agrees closely with the experimental valued 
14.26 A 3 /atom. 
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FIG. 6: Top: Energy differences AE of the fee, hep and simple hexagonal (Af) structures with 
respect to bec in Tc as function of volume. Bottom: Equation of state of Tc up to 500 GPa. 



F. Ruthenium 



At ambient pressure, the structure of Ru is hep. Experiments have been performed up 
to 56 GPa, and no phase transition has been found. Our results for the energies of the fee 
and hep structures relative to bee (Fig. [7j) show that hep is the most stable up to pressures 
of at least 400 GPa. As can be seen in Fig. [7J our calculated EOS of hep Ru is in very- 
good agreement with recent experimental data.— In particular, we obtain an equilibrium 
volume of 13.42 A 3 /atom which is close to the experimental value^ 1 ^ of 13.47 A 3 /atom. 
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FIG. 7: Top: Energy differences AE of the hep and fee structures with respect to bee in Ru 
as function of volume. Bottom: Equation of state of Ru up to P ~ 400 GPa compared with 
experimental data. 63 

(An earlier calculated equilibrium volume 6 ^ of 13.57 A 3 /atom is in rather poor agreement 
with our value.) We also note that our calculated c/a value of 1.58 is exactly the value found 
experimentally.— 
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G. Rhodium and Palladium 

The last two elements treated here, Rh and Pd, both have the fee structure at ambient 
pressure, and retain this structure up to the highest pressures reached so far experimen- 
tally^ 1 ^ (~ 50 and 77.4 GPa, respectively). Our results for the energy differences fee - bec 
and hep - bee (Figs. [8] and [9]) give no indication that any phase transition will be found in 
the range up to ~ 500 GPa. As for most of the other 4d elements, our calculated EOS results 
(Figs. [8] and [9]) are in close agreement with experiment. Our calculated equilibrium volumes 
of 13.71 and 14.73 A 3 /atom for Rh and Pd respectively agree closely with the experimental 
values*^ of 13.75 and 14.72 A 3 /atom. 

H. Zero-temperature phase diagram 

The pressures of the transitions for each element (atomic number Z) reported above can 
be regarded as points on phase boundaries drawn in the (P, Z) plane. These boundaries 
are described by the dependence of P on Z treated as a continuous, rather than a discrete 
variable. In the real world, the elements form a discrete series, and only integer values 
of Z are available, but in DFT theory there is no difficulty in treating Z as continuous. 
(In tight-binding theories of transition-metal energetics, it is common practice to treat the 
number of electrons per atom as continuous.) However, to save effort, we have not actually 
attempted to perform FP-LAPW calculations for non-integer Z, preferring to obtain the 
phase boundaries by simple interpolation. 

To perform the interpolation, we note that at T = K the enthalpies H = E + PV 
of coexisting phases must be equal. To take an example, the hep — bec energy differences 
for Mo (Z = 42) and Tc (Z = 43) at P = 400 GPa are 0.62 and -0.22 eV/atom, and 
the hep — bec volume differences are —0.14 and —0.12 A 3 /atom, so that the hep — bec 
enthalpy differences are 0.27 and —0.52 eV/atom. By linear interpolation, we estimate 
Z = 42.34 as the coexistence value between bec and hep at 400 GPa. We have used this 
interpolation scheme to estimate the bec — hep boundary between Mo and Tc, and the 
hep — fee boundary between Ru and Rh. For the complicated region Y — Zr — Nb, we have 
performed interpolation only for the fee — bec boundary passing through the Y — Zr region; 
in the other cases, the boundaries have been drawn approximately by means of straight 
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FIG. 8: Top: Energy differences AE of the fee and hep structures with respect to bee in Rh 
as function of volume. Bottom: Equation of state of Rh up to P ~ 500 GPa compared with 
experimental data,- 2 

lines. The resulting generalized phase diagram is shown in Fig. [10] • 



IV. DISCUSSION AND CONCLUSIONS 

We noted in the Introduction the general expectation that the sequence of crystal struc- 
tures across a transition-metal series will be the same at high P as at low P, but that the 
phase boundaries should slope to the left (dP/dZ < 0). The well-known ideas behind this ex- 
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FIG. 9: Top: Energy differences AE of the fee and hep structures with respect to bec in Pd 
as function of volume. Bottom: Equation of state of Pd up to P ~ 450 GPa compared with 
experimental data, 72 

pectation are that (a) the sequence of stable structures depends mainly on the band energy, 
i.e. the sum of single-electron energies; (b) the band energy depends mainly on the <i-band 
densities of states in the different crystal structures, and on the number of (i-electrons; (c) for 
given atomic number Z, increase of P causes the bottom of the <i-band to rise with respect 
to the o?-band centroid, so that the effective number of ^-electrons increases. Because of this 
last effect, increase of P should cause a given element to become more like its neighbour 
on the right. This is expected to be a stronger effect for early transition metals, since the 
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FIG. 10: Calculated zero-temperature phase diagram of the 4d transition metals series. Solid 
squares correspond to phase transitions obtained directly from the FP-LAPW calculations while 
open squares show results obtained by interpolation (see text). 

increase of the effective number of <i-electrons with compression is larger for them.— 

These expectations are confirmed most clearly by the major boundary separating bcc 
from more close-packed structures on the left of the phase diagram. This boundary causes 
Zr at rather moderate pressures to adopt the low-P structure of Nb, and Y to do the same at 
much higher pressures. However, the boundary separating bcc and hep in the middle of the 
diagram is almost vertical up to P ~ 100 GPa, and it is only at much higher pressures that 
it slopes to the left; furthermore, it appears that Mo never adopts the low-P structure of its 
neighbour Tc. Similarly, the boundary between hep and fee on the right of the diagram is 
practically vertical up to P ~ 500 GPa. 

DFT predictions of the many transitions shown by Y and Zr already give evidence that 
DFT provides reasonably good quantitative accuracy for the pressure dependence of band 
structures. The bcc-fcc transition in Mo will probably remain out of reach of experiment 
for some time. On the controversial question of whether the high-P transition in Mo is 
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really bcc-fcc or bcc-hcp, our results support the former. The reason for this controversy is 
clearly that the hep and fee energies are very close to each other, so that the calculations 
have to be converged to high accuracy to yield reliable conclusions. We have taken pains to 
ensure that our prediction of bcc-fcc is robust with respect to all convergence parameters. 
Consequently, we believe that the small and unexpected field of fee at high P and roughly 
half band filling is a real effect. 

The bcc-fcc transition in Mo is relevant to the interpretation of the shock experiments 
mentioned in the Introduction. When this transition was first studied theoretically, the 
predicted transition pressure was 420 GPa, 20 which is much lower than our value of 660 GPa. 
At that time, it seemed likely that the transition is closely related to the transition seen in 
shock experiments^ on Mo at P ~ 210 GPa and an (estimated) temperature of 3100 K. Our 
confirmation of earlier work giving a much higher transition pressure makes it much less likely 
that there is a direct connection with the shock transition. The only way of maintaining 
this connection would be to postulate that the transition pressure decreases strongly with 
increasing temperature. However, our DFT calculations of phonon frequencies in the two 
structures^ 3 - show that the transition pressure actually increases with temperature, so that 
a direct connection with the shock transition is completely ruled out. 

Turning now to detailed comparisons with experiment, we have shown that all our calcu- 
lated equations of state P(V) agree very closely with experiment. The sequences of stable 
structures with increasing pressure are always correctly predicted, where experimental data 
are available, but our calculations have a clear tendency to underestimate transition pres- 
sures Pt by typically ~ 10 GPa. This can be attributed to a relative shift 5E of the energies 
of coexisting phases. Since at T = K, the enthalpies H of coexisting phases are equal, a 
shift of Pt is related to SE by: 

5E = 5P t {{dHjdP)? - (dH 2 /dP) T ) , (1) 

where Hi(P) and H 2 (P) are the enthalpies of the two phases. Since (dH/dP)x = V, the 
shift 5P t can be estimated as SP t = SE/ (V1 — V2), where V\ and V 2 are the coexisting volumes. 
Then a shift 5Pt = 10 GPa equates to a ratio 8E/(Vi — V 2 ) = 60 meV/A 3 atom. Taking 
an example, the volume change V\ — V 2 at the hcp-a> transition of Zr was calculated to be 
0.1 A 3 /atom, so it would require a relative shift of 6 meV/atom to account for the error in 
P t . On the other hand, for the u-bec transition of Zr, for which V\ — V 2 = 0.6 A 3 /atom, 



23 



the relative shift would have to be 36 meV/atom. It is easy to see that room-temperature 
thermal effects are unlikely to be the cause. Estimating the vibrational free energy per 
atom F vi b as 3k^T ln(hu / k^T) , with uo the geometric-mean phonon frequency, a difference 
5uj between coexisting phases causes a relative free energy shift <5P v ib — ^k-^T 5uj / uj . With 
3kftT = 80 meV at room temperature, unreasonably large 5uo/uo values would be required. 
On the other hand, we have cited evidence (Table IHIj) that differences between density 
functionals can lead to P% shifts of order 10 GPa in the tu-bcc transition of Zr. However, 
there is a third possible cause, namely errors of DFT implementation. We have made efforts 
to ensure that our FP-LAPW energies are converged to ~ 1 meV/atom with respect to 
plane-wave and angular-momentum cut-offs, but there remain possible linearization errors, 
whose size is difficult to estimate. It seems clear that implementation errors must be the 
reason for some of the very large differences between DFT predictions of Pt values. We think 
this explains the difference of over 200 GPa in the predicted P t of the fcc-bcc transition in 
Y (Table [IT]) , because the earlier DFT work was done at a time when the treatment of semi- 
core states was less well developed. The range of ~ 300 GPa in predicted P t values for the 
high-P transition in Mo must also be attributed to implementation errors in early work. 

We conclude by recalling that the present calculations on the zero-temperature (P, Z) 
phase diagram are intended as a prelude to the systematic mapping of the (P, T, Z) diagram, 
including the solid-liquid coexistence surface as a function of P and Z. In spite of major 
progress in the first-principles calculation of melting curves over the past 10 years,— iZlZLZS 
the computation of the entire (P, T, Z) phase diagram is clearly a major challenge, which 
will need to be tackled in stages. The use of DFT to calculate harmonic vibration frequen- 
cies^ will allow quite rapid progress at temperatures up to about one third of the melting 
temperature, and we have already reported systematic calculations of this kind for Fe, Ta 
and Mo over a wide range of pressures.— 1 ^ For melting curves, several first-principles meth- 
ods are available.— 1 ^ However, we believe it will also be valuable to make rapid and more 
approximate surveys using tight-binding methods, and we plan to report on this approach 
in the near future. 
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